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Graphene has two atoms per unit cell with quasiparticles exhibiting the Dirac-like behavior. 
These properties lead to interband in addition to intraband optical transitions and modify the /- 
sum rule on the longitudinal conductivity. The expected dependence of the corresponding spectral 
weight on the applied gate voltage V g in a field effect graphene transistor is ~ const — |Vg| 3//2 . For 
V g — 0, its temperature dependence is T 3 rather than the usual T 2 . For the Hall conductivity, the 
corresponding spectral weight is determined by the Hall frequency ljh which is linear in the carrier 
imbalance density p, and hence proportional to V g , and is different from the cyclotron frequency for 
Dirac quasiparticles. 

PACS numbers: 78.20.Ls, 74.25.Gz, 73.50.-h, 81.05.Uw 



I. INTRODUCTION 

The real part of frequency dependent optical conductivity a xx (Q) is its absorptive part and its spectral weight 
distribution as a function of energy (Ml) is encoded with information on the nature of the possible electronic tran- 
sitions resulting from the absorption of a photon. Even though the relationship of the conductivity to the electronic 
structure and transport lifetimes is not straightforward, much valuable information can be obtained from such data. 
In particular, the /-sum rule on the real part of a xx {VL) stated in its simplest form for an infinite free electron band, 



<KlRea xx (Q.) = -f (1.1) 
o 

has proved particularly useful. Here top is the plasma frequency, up = 4irne 2 /m, with n the free carrier density 
per unit volume, — e < the charge of electron, and m its effective mass. In this case the right hand side (RHS) of 
Eq. (|1.1[) is independent of temperature and of the interactions. More generally for a finite tight binding band, the 
optical sum rule takes the form 

Qm e 2 v-^ <9 2 ek 

dflRea xx (Cl) = — yn^ a —— , (1.2) 
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with Qm a cutoff energy on the band of interest and only the contribution to Reo- xx (Q) of this particular band is to 
be included in the integral. Here h is Planck's constant, V is the crystal volume, a is the spin, ek is the electronic 
dispersion, k is the wave vector in the Brillouin zone, and nk. CT is the probability of occupation of the state |k, a). 
If we assume a free electron dispersion, et = ft 2 k 2 /2m, Eq. (|1.2[) immediately reduces to Eq. (|1.1[) with n = N/V, 
where N = ^ k nk,<r is the total number of electrons in the band. For tight-binding dispersion with nearest neighbor 
hopping on a square lattice it is easy to show that the RHS of Eq. (| 1 .2(1 reduces to e 2 /h 2 multiplied by minus one-half 
of the kinetic energy, Wk.e. per atom. In this particular case the optical sum rule can be used to probe the change 
in kinetic energy of the electrons as a function of temperature or with opening of a gap in the spectrum of the 
quasiparticle excitations which has been a topic of much recent research [l|, [H, 0, [B[ . 

When a constant external magnetic field B is applied to a metallic system in the z-direction, the optical conductivity 
acquires a transverse component a xy (Q) in addition to the longitudinal component a xx (Q). This quantity gives 
additional information on the electronic properties modified by the magnetic field. In this case Drew and Coleman 
Q have derived a new sum rule on the optical Hall angle 0h(Q). If we define 

t H (n) = tzn9 H (n) = ?^, (1.3) 
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then 
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(1.4) 



where ojh is the Hall frequency which corresponds to the cyclotron frequency u> c = eB/mc for free electrons. Here c 
is the velocity of light. 

Recently, graphene which is a single atomic layer of graphite, has been isolated Q and studied 

HQ (seeRefs.OjJEI 

for a review). The two-dimensional (2D) graphene honeycomb lattice has two atoms per unit cell. Its tight binding 
band structure consists of the conduction and valence tt bands. These two bands touch each other and cross the Fermi 
level, corresponding to zero chemical potential [i, in six K points located at the corners of the hexagonal 2D Brillouin 
zone, but only two of them are inequivalcnt. The extended rhombic Brillouin zone can be chosen to contain only 
these two points inside the zone. In a field effect graphene device 0, H, [1] electrons can be introduced in the empty 
conduction band or holes in the filled valence band through the application of a gate voltage and thus the normally 
zero value of the chemical potential fi can be changed continuously. In this paper we wish to consider the optical 
sum rules described above for the specific case of graphene. Two essential modifications arise. The existence of the 
two bands means that interband transitions need to be accounted for in addition to the usual intraband transitions 
of the previous discussion. Second, the electronic dispersion curves near K points are linear in momentum and the 
low-energy quasiparticle excitations are described by the "relativistic" (2 + 1) dimensional Dirac theory [l2l. [l3l [bil [l5| 
rather than the Schrddinger equation. 

The effective low-energy Dirac description of graphene turns out to be insufficient for the derivation of the sum 
rules. This can be easily understood from two examples. The RHS of Eq. (| 1 . 2[) is normally related (Refs.0,0, H, [H) to 
the diamagnetic term that is defined as the second derivative of the Hamiltonian [see Eq. (|2.7p below] with respect to 
the vector potential. This term is zero within the Dirac approximation. On the other hand, it follows from the same 
Dirac approximation that, in the high-frequency limit, the interband contribution to conductivity is constant [l6l. [l7j 



The frequency f2 is unbound in the Dirac approximation, although physically f2 should be well below the band edge. 
This example indicates that when considering sum rules one should go beyond the Dirac approximation. 

The second example is that the cyclotron frequency oj c for the Dirac quasiparticles [19] is defined in a different way, 
oj c = eBv%/Jc\iJ,\), where Vf is the Fermi velocity. This definition follows from the fact that a fictitious relativistic 
mass [H 0, US] m c = \uU v% p lays the role of the cyclotron mass in the temperature factor of the Lifshits-Kosevich 
formula for graphene [2ll . |22| . This cyclotron frequency diverges as fi —* 0, also posing the question as what one 
should use as a Hall frequency on the RHS of Eq. (|1.4[) . 

It turns out that these problems are resolved when a tight-binding model is considered, but still due to the specifics 
of graphene, the /-sum rule cannot be guessed from Eq. (|1.2|) applied for the two-band case. Further important 
modifications to the /-sum and Hall angle sum rule for graphene are expected and found in this paper. 

The paper is organized as follows. In Sec. [IT] we introduce the tight-binding model and the necessary formalism. 
In Sees. IIIII and ITVl we describe the derivation of the sum rules for the diagonal and Hall conductivities and consider 
the specific case of graphene. The partial sum rule for the Hall angle is investigated in Sec. [V] In Sec. IVI1 the main 
results of the paper are summarized. Mathematical details are given in two Appendices. 

II. TIGHT-BINDING MODEL AND GENERAL REPRESENTATION FOR ELECTRICAL 

CONDUCTIVITY 



The honeycomb lattice can be described in terms of two triangular sublattices, A and B. The unit vectors of the 
underlying triangular sublattice are chosen to be 



(see also Ref. 



Rea xx (fl) ~ — , 



n > m, t. 



(1.5) 



A. Tight-binding model: Paramagnetic and diamagnetic parts 



ai = (aV3/2,a/2), 
a 2 = (aV3/2, -a/2) 



(2.1) 



where the lattice constant a = |ai| = \a.a\ = \/3acc and acc is the distance between two nearest carbon atoms. Any 
A atom at the position n = aini + a^ni-, where n\,ni are integers, is connected to its nearest neighbors on B sites 
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by the three vectors 6i 



5i = -(ai + a 2 )/3, 
6 2 = 2/3ai - l/3a 2 , 
5 3 = -l/3ai+2/3a 2 . 

We start with the simplest tight-binding description for 7r orbitals of carbon in terms of the Hamiltonian 



a n,<r ex P ( fc^ A ) & n+<S,<r + C.C. 



(2.2) 



(2.3) 



where t is the hopping parameter, a n/T and b n+ s,<r are the Fermi operators of electrons with spin a on A and i? 
sublattices, respectively. Since we are interested in the current response, the vector potential A is introduced in the 
Hamiltonian (|2.3|) by means of the Peierls substitution aj, a b m ,a ->«J, ff exp (—It / Adr) 6 miCr , which introduces the 
phase factor exp(||<5A) in the hopping term (see Ref. [1 for a review). We keep the Planck constant h and the velocity 
of light c, but set ks = 1. 

Expanding the Hamiltonian (|2.3[) to the second order in the vector potential, one has 



-A(n)j(n) - ^^(^^(11)^(11) 



a, /3 = 1,2. 



The total current density operator is obtained by differentiating the last equation with respect to A a (n), 

dH 



Ja(n) 



and consists of the usual paramagnetic part 



d(A a /c) 



la ( n ) = ~Y^2 Sa [ a kv b "+S,<r ~ b l+S,a a n,, 



j a (n)-r afi (n)Ap(n)/c, 



(2.4) 



(2.5) 



(2.6) 



and diamagnetic part 



T Q( 3(n) 



d 2 H 



te z 



d(A a /c)d(A /c) V 



S.a 



(2.7) 



B. Noninteracting Hamiltonian 

The noninteracting Hamiltonian Hq in Eq. ()2.4[) written in the momentum representation reads 



H " = J2[ 7f^Tt(k)H T CT (k), Ho 



with 0(k) =tJ2g e ikS = e(k)e lv ^ . The dispersion law e(k) is 



e(k)e 4v(k ) 
e(k)e~^ k > 



e(k) = J 1 + 4 cos cos ^ + 4cos2 



In Eq. (|2.8[) we introduced the spinors 



T.(k) 



a CT (k) 
6 CT (k) 



with the operator T cr (k) being the Fourier transform of the spinor T CT (n) = y ^ 

T tT (n) = VS [ -^e ikn T CT (k) 

JBZ { Z7r ) 



^n.cr 



(2. 



(2.9) 



(2.10) 



(2.11) 
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Here S — v / 3a 2 /2 is the area of a unit cell and the integration in Eqs. (|2.8[) and (|2.11[) goes over the extended 
rhombic Brillouin zone (BZ) which is characterized by the reciprocal lattice vectors t>i = 2ir /a(l/ y/3, 1) and b2 = 
27r/o(l/V5, -1). 
We also add the term 

Ho^Ho-fiY'f -^-Tt(k)/T CT (k) (2.12) 



BZ 



(27T) 



with the chemical potential /i to the Hamiltonian H , so that our subsequent consideration is based on the grand 
canonical ensemble. The corresponding imaginary time Green's function (GF) is defined as a thermal average 

G4n-T 2 ,n 1 -n 2 ) = -(T r T .(n,ni)T+(75,iia)) (2-13) 

and its Fourier transform is 

Gain -r 2 ,ni - n 2 ) =TV / . G(w K , k) exp[-w„(ri - t 2 ) + ik(ni - n 2 )], (2-14) 

with 

nc \*\ ( iuJn + / i K + cr + ( M k ) +<T_0*(k) 

G(iw n ,k) = ■ — -2 ^-r , w n = 7r(2n + 1)T, (2.15) 

where the matrix <r± = [a\ ±za 2 )/2 made from Pauli matrices operates in the sublattice space. In Eq. (|2.15[) the spin 
label <7 is omitted, because in what follows we neglect the Zeeman splitting and include a factor of 2 when necessary. 
The GF (|2.15[) describes the electronlike and holelike excitations with energies £'±(k) = ±e(k) — fi, respectively. The 
dispersion e(k) near K points is linear, E±(p) = ±hvF\Jp\ +v\~ Mi where the wave vector p = (pi,p 2 ) is now 
measured from the K points and the Fermi velocity is vp — \/3ta/(2h). Its experimental value [1, @| is vp w 10 6 m/s. 



C. Electrical conductivity 

The frequency-dependent electrical conductivity tensor a a p (f2) is calculated using the Kubo formula 0, Ej[ 

°<*?W- - i{Sl + i0) > K aP (Sl + tQ)=—+ — , (2.16) 

where the retarded correlation function for currents is given by 

/CO 
dte int U^(t), H« (t)=-i8(t)Tr(p[J a (t) 1 J fl {0)]), t 2 - 1 ?) 
-oo 

V is the volume of the system, p = exp(— /3-f?o) JZ is the density matrix of the grand canonical ensemble, = 1/T 
is the inverse temperature, Z = Tr exp(— (3Hq) is the partition function, and J a are the total paramagnetic current 
operators with 

J a (t) = e im / h J a (0)e- im / h , J a (t) = £ j£(t,n), (2.18) 

expressed via the paramagnetic current density (|2.6p . Using the representation for Tl a p(Q) in terms of the matrix 
elements of the current operator J a (t = 0), one can find [231 ] the high-frequency, f2 — > oo asymptotic of the Hall 
conductivity, 

**yV) = {([J„ Jy}) + j^([[{J*,H},H}J y ]) + O (^_) } , (2.19) 

while the longitudinal conductivity is 

1 f I 



Here [,] is the commutator. Below in Sees. Ill and IV we will use Eqs. (2.16), (2.19). and (2.20) to outline the formal 
derivation of the sum rules. 
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III. DIAGONAL OPTICAL CONDUCTIVITY SUM RULE 



The optical conductivity sum rule is a consequence of gauge invariance and causality. Gauge invariance dictates 
the way that the vector potential enters Eq. (|2.3p and, respectively, determines the diamagnetic and paramagnetic 
terms in the expansion (|2.4[) as well as the form of Kubo formula (I2.16|) . The causality implies that the conductivity, 
Eq. (|2. 16[) satisfies the Kramers-Kronig (KK) relation 

m J_ QO w - 11 

Combining together the high-frequency limit of the KK relation (|3.1[) . 

1 f°° 

Ima a a(n) = — P / dujRea al3 (u)), -> oo, (3.2) 

and the asymptotic (|2.20p . we arrive at the sum rule 

oo 

- [ dfiRecr^) = %l (3.3) 

TT J V 

— oo 

Taking into account that Rea xx (£l) is an even function of O, we observe that for a single tight-binding band {T a p)/V 
corresponds to the RHS of Eq. (|1.2[) . Also defining the plasma frequency ujp via ujp/(4:ir)8 a f3 = (Tap) /V we can 
rewrite the sum rule (|3 . 3[) in the form Below we calculate this term for the Hamiltonian (|2.3|) with a nearest 

neighbor hopping on the hexagonal graphene lattice. 

A. Explicit form of the diamagnetic term 

The diamagnetic or stress tensor (r Q| a) in the Kubo formula (|2.16[) is a thermal average of Eq. (|2.T[) 

(rctf) = <$>«/?(*)>■ (3-4) 

n 

This term is calculated in Appendix [A] and is given by 

e(k), (3.5) 



2<= 2 f d 2 k 

" ' W(e(k))-n F (-e(k))] 



V ft 2 / HZ (2tt) 2 



dfc 2 ^ dk a ) 



where np(uj) — l/[exp((u> — fJ-)/T) + 1] is the Fermi distribution. Because the graphene structure contains two atoms 
per unit cell (two sublattices), there are two bands in the BZ which correspond to positive- and negative-energy Dirac 
cones. The momentum integration in Eq. (|3.5[) is over the entire BZ and the thermal factors njr(e(k)) and np{— e(k)) 
refer to the upper and lower Dirac cones, respectively. We note that a simple generalization of Eq. (|1.2|) for a two 
band case would miss the term with the derivative of the phase, (dip(k)/dk a ) . This term occurred due to the fact 
that the Peierls substitution was made in the initial Hamiltonians (|2.3p and (|2.8|) rather than after the diagonalization 
of Eq. (|2.8|) . The second comment on Eq. (|3.5[) [see also Eq. (|A4p in Appendix |A] is that (r aa ) vanishes if e(k) is 
taken in the linear approximation. This reflects the absence of the diamagnetic term in the Dirac approximation. The 
correct way is firstly to take the derivatives in Eq. (|3.5[) . This is done in Eq. (|A5|) in Appendix [X] and leads to the 
final result 



d 2 k 

V ~ 3ft 2 J BZ 



[n F (e(k))-r^(-e(k))]e(k). (3.6) 



Eq. (|3.6[) is equivalent to Eq. (|3.5[) . Note that (r QCe ) is always positive and does not depend on the arbitrary choice 
of the sign before t in Eq. (|2.3p . Now Eq. (|3.6p is e 2 /h 2 times — 2/(3\/3)(^ —0.39) of the kinetic energy per atom 
instead of —1/2 for the usual square lattice. At zero temperature for /x > (to be specific) the lower Dirac cone is 
full and, in the conductivity, only interband transitions are possible for these electrons. Similarly, the electrons in 
the upper Dirac cone can undergo only intraband transitions. In the above sense the first thermal factor in Eq. 
corresponds to intraband and the second to interband at T = 0. 
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It is useful to separate explicitly the contribution (T xx (fi = T = 0)) of the Dirac sea from Eq. (|3.6p : 

(t xx ) = (t xx (jjl = T = 0)) + (t*(jjl,T)), (3.7) 



where 

(r^(/i = T 1 = 0)) e 2 a 2 /" d 2 k 



1/ 3ft 2 7 SZ (2tt) 2 

is the contribution of the Dirac sea (the energy of the filled valence band) and 

(t£(h,T)) eV f d 2 k 



(-e(k)) (3.8) 



V 3fr 2 



/ 7y4 M<k)) + 1 - "F(-e(k))] e(k) (3.9) 
Jbz { Z7T j 



is the electron-hole contribution. Expressions (13. 6p . (13. 8|l . and Ij3.9|) contain the energy e(k) as happens also for the 
case of the square lattice with nearest neighbor hopping mentioned in the Introduction. 

The numerical calculation of the Dirac sea contribution (|3.8I) with the full dispersion (|2.9p gives 

(t xx (vl = T = 0)) eH 

— = a— T , a « U.ol. (3.10) 

V nr 

The same answer also follows from the linearized Dirac approximation with the trigonal density of states if the band 
width W is taken to be W = \J \/3nt. The electron-hole contribution (|3.9|) can be estimated analytically in the linear 
approximation for the dispersion law, 



V Znh^vl 



Li 3 (-e"/ T ) +Li 3 (-e-^ T )l , (3.11) 



where Lis(z) is the polylogarithmic function [24J. This shows that for n = the temperature dependence of the 
diagonal conductivity sum rule is ~ T 3 , in contrast to a well-known T 2 dependence 0, El, 0, HI- We note, however, 
that because (t^((i — 0,T))/V ~ —T 3 /t 2 which is small, the T 3 behavior is unlikely to be observed. On the other 
hand, using asymptotics of the polylogarithmic function for \fi\ 3> T, Eq. p. lip can be written in the form 



V 97rn 4 w 2 



F 



[H 3 + tt 2 HT 2 ]. (3.12) 



The |/x| 3 behavior at T — is easily understood from Eq. (|3.9p in which case it reduces to an integral over energy 
ranging from to /i of the density of states (DOS) which is proportional to |e| in graphene. The integrand is therefore 
proportional to e 2 leading directly to the |/x| 3 dependence. This cubic power law reflects directly the Dirac nature 
of the electronic dispersion relation encoded in the linear in the e DOS. At finite T the DOS in Eq. (|3.9p for fi = 
provides an additional factor of T as compared with the conventional case, leading to T 3 rather than T 2 behavior. On 
the other hand, for \/j,\ ^> T the DOS can be evaluated at fi, leading to a |/i|T 2 law. All these results follow directly 
from a linear dispersion of the massless Dirac quasiparticles. 

We note that the |/x| 3 behavior is more likely to be observed than the T 3 dependence by varying the gate voltage 
V g . Using Eq. (|4.7|) for carrier imbalance p which is proportional to V g , we obtain that (t xx (V 9 )) /V ~ — |K,| 3 / 2 . Note 
that as one can see from Eq. (|3.12p . the electron-hole contribution (r^fp, T)) in Eq. (|3.9p is negative and it has to 
be added to Eq. (|3.8p to get the total non-negative contribution to the optical sum. If at T — we take the chemical 
potential to fall at the top of the band, so that the positive energy Dirac cone is completely occupied, Eq. (|3.9p reduces 
to Eq. (|3.8p except for a difference in sign. The two contributions cancel, because a full band cannot absorb. The 

same result follows from Eq. (|3.12p when [i is taken to be the energy cutoff on the band-width W = \J y/Sirt. We also 
note from Eq. (|3.12p that for 3> T we recover the already mentioned T 2 law [1, [1, H H . 

IV. HALL ANGLE SUM RULE 

We now consider the optical Hall angle (|1.3p that plays the role of the response function to an injected current 
rather than an applied field: 

j x (il) = a xy (n)E y (n) = t H (Q)j y (Cl), j y (Q) = a xx {Cl)E y (Q). (4.1) 
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It was proved in Ref. @ that the response function tjj(fi) satisfies the KK relation 

t H (Sl) = — P / duj—^—t H {uj). (4.2) 
m J_ oc LU-tt 



Multiplying Eq. (14. 2p by — ifl and taking the limit £1 — > oo, we obtain the sum rule 

1 

with the Hall frequency 



du>t h (to) — oj h , (4-3) 



cjh = lim [-iOte(O)]- (4.4) 
|n|— >oo 

The Hall angle sum rule (|1.4|) follows from Eq. (|4.3p after we take into account that the real and imaginary parts of 
tni'jj) are even and odd functions of f2, respectively. Microscopic considerations based on the Kubo formula (|2.16|) 
show that the high-frequency limit (|4.4[) exists and is given by 

lim ^v(n) = ([^ = o),y = Q)]) > (45) 

where in the last equality we used the high-frequency asymptotics, Eqs. (|2.19|) and (|2.20p . 
The commutator I XtV = ([J x (t = 0), J y (t = 0)]) is calculated in Appendix IB] where we obtain 

e 2 a 4 t 2 heB . , 

' y = ~ l — i¥c — exyVp - ( 6 ) 

Here p is the carrier imbalance (p — n e — n^, where n e and rih are the densities of electrons and holes, respectively), 
and e a b is antisymmetric tensor. The carrier imbalance for B = T = and in the absence of impurities is 

P = 1M£- (47) 

Substituting Eqs. (14^) and (|XTU)l in (|4.5p we finally obtain 

1 e5ta 2 , 

Since to? /h 2 has the dimensionality of the inverse mass and pa 2 is dimensionless, Eq. (|4.8p has the correct dimension- 
ality of the cyclotron frequency. Substituting Eq. (|4.7p and expressing vp via i, one can rewrite 

^ = -^f^ 2 (B)^. (4.9) 
97ra /it 



Here L(B) = \ eB \ hv 2 F / c is the Landau scale which in temperature units is equal to L 2 (B)[K 2 ] — 8.85 x 
10" 8 K 2 w|,(m/s)B(r) and h = 7.638 x IQ- 12 K ■ s. 

As mentioned in the Introduction, in the recent interpretation of Shubnikov de Haas measurements a gate voltage- 
dependent cyclotron mass was introduced [1, [§] through the relationship \p,\ = m c vp. If this is used in Eq. (|4.9p we 
get 

eB f p \3 



-—(^-). (4.10) 
cm c \ l.62tJ y ' 



Since a full upper Dirac band corresponds to a value p — W — y y/Sirt ~ 2.33* (see the end of Sec. IIII Ap . in this 
case formula (|4. 1 0[) resembles the formula uh = u c = eB/mc from Ref. 6 for two-dimensional electron gas with m c 
replacing the free electron mass. In graphene, however, m c varies as the square root of the carrier imbalance \p\ and 
the two cases look the same only formally. 

Finally, we notice that when the spectrum becomes gapped with E — \ffi 2 v 2 F (p\ + p?,) + A 2 , the carrier imbalance 

is 

P= = ^(M 2 -A 2 M M 2 -A 2 )sgn / . (4.11) 

V F 



This implies that the gap A can be extracted from the change in ujh obtained from magneto-optical measurements. 
This kind of measurement which reveals gapped behavior has already been done on the underdoped high-temperature 
superconductor YBa^CusOg+a; [25j. The recent measurements done in epitaxial graphite [26] and in highly oriented 
pyrolytic graphite [27J lead us to expect that this experiment should be possible for graphene. Below we restrict 
ourselves to the A = case. 
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V. PARTIAL SPECTRAL WEIGHT FOR THE HALL ANGLE SUM 



So far we have considered only the complete optical sum involving integration of the optical spectral weight over 
the entire band. The temperature dependence of such a sum has been central to many recent studies [!, 0, Q focused 
on the possibility of kinetic energy driven superconductivity in the cuprates. The effects are small and the experiment 
is difficult. The conductivity is needed up to high energies and there is no unambiguous criterion to decide where 
the band of interest may end. On the other hand the partial optical sum to a definite upper limit Sl m can give 
important information on the approach to the complete optical sum. It has also proved very useful in understanding 
the spectral weight redistribution with temperature or phase transition (see. e.g., the recent works in Refs. 28 and 
l29h . A discussion of the effect of an external magnetic field on the partial optical sum (or weight of the main peak) 
for the longitudinal conductivity in epitaxial graphite can be found in Ref. |26j. While the optical Hall conductivity is 
not as widely studied as the longitudinal conductivity, many new studies (see, e.g., Ref. have shown its usefulness 
as it gives information on the change in microscopic interactions brought about by B. The available experimental 
data for graphite [27| already contain information about optical Hall conductivity, and when the same measurements 
are done on graphene, they can be directly compared with the results discussed below. 

We saw in our previous paper [TtJ that the longitudinal conductivity 

W xx (fl m ) = / dnRea xx (n) (5.1) 
Jo 

showed (see Fig. 7 in Ref.[l7t) plateaux with the steps corresponding to the various peaks in the diagonal conductivity. 
Here we are interested in the corresponding quantity associated with the Hall angle Reijj(O) and consider the weight 

W(Cl m )= dilRet H (n). (5.2) 
Jo 

It is instructive to begin with a discussion of the frequency dependence of Reiff(Sl). As already mentioned after 
Eq. (|4.3p . because the imaginary part of both a xy (Tt) and cr xx (f2) is an odd function of SI, only the real part of 
t// (SI) contributes to the sum rule (jl.4p . Thus we will need to consider only Rei#(Sl) defined by Eq. (|1.3p . General 
expressions for the complex conductivities a xx (£l) and a xy (Tl), derived in the Dirac approximation, are given by 
Eqs. (9), (11) and (10), (12), respectively, of our previous paper [l7| . 

In what follows we will consider explicitly possible experimental configurations. For a fixed cutoff (Sl m ) sweeping 
the magnetic field gives direct information on the change in optical spectral weight, in this frequency region, brought 
about by B. For the field effect transistor configuration used in Refs. 0, H, the chemical potential fj, can be changed 
through an adjustment of the gate voltage at fixed B and S7 m . Another configuration that might be considered is to 
fix B and \i and change Sl m . These three cases serve to illustrate what can be expected in experiments. 

In Fig. [1] we show results for the frequency dependence of Re i# (SI) for a case with external magnetic field B = 0.4 T, 
at temperature T — 5K and impurity scattering rate T = 10 K. Accordingly, the energies of the Landau levels are 
M n = y/2nL 2 (B) [see also Eq. (5) in Ref. [T^— viz., M x ~ 265 K, M 2 ~ 375 K, M 3 ~ 460 K, and M 4 ~ 531 K, 
respectively. Thus the long dashed (red) curve with /i = — 20 K is for \fi\ < Mi, the dash-dotted (black) curve with 
fi = -300 K is for Ml < < M 2 , the solid (blue) line with ^ = -400 K is for M 2 < |/u| < M3 and the short dashed 
(green) curve with /i = — 500 K is for M3 < < M4. The values chosen for all the parameters used in Fig. Q]are 
quite reasonable. For a magnetic field 0.4T the various Landau energies fall within the range of available optical 
spectrometers and this is one value of B used in Ref. [26| on the optics of several-layer graphite. The value of the 
broadening parameter T is not so well known, but our choice is reasonable when compared with the line broadening 
seen in experiments. [26| Higher values would simply broaden the lines seen in Fig. [I] Further in field-effect graphene 
devices 0, B the chemical potential can be easily varied from — 3600 K to 3600 K simply by changing the gate 
voltage. Low temperatures can also be easily achieved. Starting with the long dashed (red) curve of Fig. Q] we 
notice a Drude-like behavior of Rei#(Sl) at small SI, followed by a region where it is small. It then changes sign, 
after which it exhibits a negative peak around 375 cm -1 . This is followed by a further series of small peaks, all 
negative, which decay in amplitude as the frequency S7 increases. The frequencies of the peaks do not correspond to 
the sum or differences of two adjacent Landau level energies M n as do the peaks in both Hea xx (Q) and Im <r xy (£l) 
seen in Ref. [l?]. Rather they fall, approximately, at energies between two peaks in Rea xx (fl) (see Fig. 1 in Ref. [13 
taking into account that it is plotted for B = IT), where the Im a xx (SI) also crosses zero and the denominator in 
Hetji — Re(a xy / 'o xx ) which is equal to (Re a xx (S7)) 2 + (Ima xx (fl)) 2 consequently has a minimum. The numerator 
a xx (Q)Ke a xy (fl) + Ima xx (fl)Ima xy (fl) plays little role in the peak positions, but determines their sign which can 
be positive or negative as seen in the other three curves (for larger values of chemical potential). In this case the first 
peak is positive and this is followed by, decaying in amplitude, a sequence of higher-energy peaks, all of which are 
negative. Note that the position of the first positive peak in each of the dash-dotted, solid, and dashed curves moves 
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progressively to higher energy as the chemical potential crosses a new Landau level. This happens because, as we have 
already stated, the peaks in Rei#(f2) fall in the between two consecutive peaks in Re a xx (Q) (for Mat < \fi\ < M/v+i, 
the first interband peak in Rea xx (£l) falls at Mm + M^+i)- 

Returning to the small- f2 behavior of Ret#(0), where a Drude-like behavior was noted, we see that the weight 
under this peak is progressively depleted as /i crosses through larger Landau level energies. These features can be 
understood qualitatively from approximate expressions for Rei# (f2) obtained by keeping only the leading term in the 
formulas for both the longitudinal a xx (Q) and Hall cr xy (Q) conductivity. In Ref . [l7l we provided expressions (11) and 
(12) which we will not repeat here. Both involve a sum over Landau levels n = 0, 1, 2, . . .. If \n\ falls in the interval 
]Mq, M±[, retaining only the n = contribution to the sum, we get for T = 



e 2 v 2 F \eB\ 2{Sl + 2iY) 

Mi Ml -(n + 2iT) 2 



and 



e 2 v 2 7 eB sgn /i 

~Vc M 2 - (O + 2iT) 2 ' 



^j/(") = iL ~Z: 772 /-r. , n.-F^ ' ( 5 - 4 ) 



It follows from Eqs. (pT4")l and (fOj) that 



iy / 2geg|g7csgn (eg) sgn/x Mi2r S gn(e5) sgn^ 

This shows that while Rea xx (£l) and Im a xy (Vt) are peaked near Q, = Mi, the small f2 limit of Ret^(^) is indeed 
Drude in shape with width 2L. Also its height at VI = and fi < is M±/2T which is set by the inverse scattering 
rate 2r and by the Landau scale L(B), because Mi = \/2L(B). 

On the other hand, for /i g]Mjv,Mjv+i[ the first nonzero term in the expression for the Hall angle tji(Cl) is for 
n = N . Retaining only this contribution to the conductivities (11) and (12) from Ref. 17 and defining 

A " = Jmn^Tm^y^Js^^Y (5 ' 6) 



for T = we get 



_ e 2 v 2 F \eB\{n + 2iT) ( A N A+ \ 

M ' ~ tt« {m n+1 - M N + M N+1 + M N J ' ^ 

a xy (Q) = -i^l( A J f +A+). (5.8) 
Hence for tu we obtain the Drude term 

^sgn(eg)w A N + A+ 

n + 2iT A N /(M N+1 -M N )+A+/{M N+1 + M N ) 

It follows directly from this formula that Hetn in the limit of large N is given approximately by 

1 Mi2T sgn (e.B) sgn tt 

Re * g( ^-2Vf n'+V ■ iV>>L (5 - 10) 



The value Ret/f(ri = 0) = Mi/(4y/NT) is roughly verified in our numerical work, where larger N corresponds to 
larger chemical potential, and these results serve as a guide to our numerical work. 

Our theory predicts not only the position in energy of the various lines as well as the shape of the Rei#(f2) at 
small photon energy, but also provides values for the optical spectral weight under the various features seen in Fig. [1] 
Information on this spectral weight is conveniently presented in terms of the partial sum rule of Eq. (|5.2p for typical 
values of the cutoff fl m as a function of the magnetic field. In Fig. [2] we show numerical results for the Hall angle 
spectral weight W{Vt m ) in cm -1 as a function of the value of the external magnetic field B for four values of the 
cutoff Q, m at fixed value of /i = — 20 K, T = IK and T = 15 K. For the long dashed (red) curve with the cutoff 
f2 m = 300 cm -1 only the frequency region in the long dashed (red) curve of Fig. [T] which falls below the first peak 
(negative in this case) in Re£#(f2) is integrated in the weight. While in Fig. [T] the field B = 0.4 T, the peaks in 
Ret#(f2) move to higher energies as B increases and so for all values of B used in the figure only the Drude region 
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around small 17 of the curve is integrated over. Using the formula (|5.5p as a rough approximation we expect in this 
case W(f2 m ) to scale as the square root of B coming from the M 1 factor. This is verified to good accuracy in the 
numerical work. When Q, m increases, as is the case for the other three curves of Fig. [2j the peaks in Rei# (fi), which 
are negative for /i = — 20 K, start entering the integral and this reduces the value of the optical sum. However, we 
have no simple approximate analytic formula which might capture the essence of the situation in this case, so we 
must rely on the numerical work. It is also clear that for a given value of B, the reduction in the optical sum caused 
by the negative peaks in Ret#(fi) will be less as B increases, because the peaks move to higher energies. In fact the 
upward steps, seen most clearly in the dash-dotted (black) curve of Fig. [5]with the cutoff f2 OT = 700 cm -1 , correspond 
to values of B for which a peak in Rei#(fi) is starting to fall outside the integration range. The first step occurs 
around B ~ 0.75 T when the second negative peak in the long dashed curve of Fig. Q]is moving through 700 cm -1 , 
while for B ~ 1.6 T it is the first peak in Rei#(fi) which is involved. This peak is larger in absolute value and so 
this second step is larger. Note also that in this case the dash-dotted (black) curve merges with the long dashed (red) 
curve for the smaller cutoff as it must, since in both instances only the Drude-like contribution at small Q is relevant 
to the integral. 

Another feature of the curves in Fig. [2] is worth commenting on. We note that the short dashed (green) curve for 
fl m — 1200 cm -1 becomes linear in B at small B in contrast to the long dashed (red) curve which, as we saw, went like 
y/B. We can get some understanding of this crossover which corresponds to the case when many peaks are involved 
in the integration, by considering the B — > limit. In Ref. IT?] we also obtained the expressions (13) for a xx and (14) 
for a xy which for A = read 



2ie 2 {fl + 2iT) 
h 



1 



(n + 2iry 



duj\uj\ 



f dnp(u 



\ duo 



°° n F (-uj) -n F (uj) 

WjJ —— — — 

o (Q + 2iT) 2 - Au 2 



dni 



(5.11) 



and 



a xy (fl) 



e 2 v F eB 



For T = we get 



( dnp{uj) dnp(-uj) 



7rc J \ dui 
1 



8u> 



(fi + 2iT) 2 Alu 2 - (fi + 2iT) 2 



(5.12) 



2ie 2 
h 



n + 2ir 



? eB 



TIC 



SgB.fl 



i 2H-(n- 
4 n 2\fi\ + (n- 
i 



2iT) 



2iT) 
1 



_(n + 2iT) 2 (n + 2ir) 2 - a^ 2 



(5.13) 
(5.14) 



Hence for ^> fi, T we obtain 



1 



heB sgn [i 



1 



Cfl 



2T-iQ 



2r-in' 



(5.15) 



where in the second equality we introduced the cyclotron mass m c . This form is to be contrasted to that obtained in 
Eq. (|5.5[) for large magnetic field B, where tjy (SI) is proportional to \f~B rather than to B as in Eq. (|5.15[) . Except 
for a sign change this result is of the same form as in Eq. (22) of Ref. @ with m c playing the role of mass which in 
graphene varies as the square root of the carrier imbalance |p|@, 1|. Here we emphasize the linear dependence on 
magnetic field B in Eq. (|5. 15|) as well as in the full sum rule (II. 4|) with loh given by Eq. (|4.9p . 

An experimental configuration which has already been used by Li et al. in Ref. |3Gj is incorporating the specimen 
into a field-effect deviceQ, 0, 0. In this case the chemical potential /j, is easily changed by varying the gate voltage. 
In Fig. [3] we show numerical results for W(Ct m ) in cm -1 as a function of chemical potential /i in K for four values of 
the cutoff f2 m . The magnetic field B has been set to 0.4 T, T = 1 K, and T — 15 K. We see steps occurring in these 
curves as the chemical potential fi crosses the Landau level energies - viz., Mi ~ 265 K, M 2 ~ 375 K, Af 3 ~ 460 K 
and M4 ~ 531 K, respectively. Between successive sharp rises, W(O m ) stays nearly constant quite independent of fi. 
The long dashed (red) curve has the smallest cutoff equal to 300 cm -1 which falls below the first peak in the long 
dashed (red) curve of Fig. [1] for \i = — 20 K and its value ~ 135 cm -1 corresponds to the area under the Drude part 
of the curve (in Fig.[T]). Remaining with Q, m = 300 cm -1 as \n\ is increased beyond M\ but is less than M2, the value 
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of W(Cl m ) decreases. This now corresponds to the dash-dotted (black) curve of Fig. Q] which has a smaller Drude-like 
piece than does the long dashed (red) curve and the drop is close to a factor of 0.56 while our simplified but analytic 
formula (|5.10p predicts 0.71 (which is better than can be expected as it is valid only for N ^> 1 and we are using it 
outside the range of validity). The dash-dotted (black) curve has fi m = 450 cm -1 which is chosen to fall above the 
first (negative) peak in the long dashed (red) curve of Fig. [TJ For \/i\ < Mi the dash-dotted curve of Fig. Q] applies, 
and to get W(f2 m ), we are integrating through the first negative peak in addition to the Drude-like peak centered 
at = 0. This reduces its value to about half the value it had for the lower cutoff fi TO = 300 cm -1 . In the region 
/j, G]Afi, Af2[, however, the dash-dotted curve of Fig. [T] applies and we are integrating over the first positive peak in 
Rei#(r2) as well as through the Drude-like region and so the value of the weight W(il m ) is now increased. After this 
it decreases again as more negative peaks are integrated over. The other two curves can be understood as well with 
arguments similar to those used for the first two cutoffs. 

Finally in Fig. Q] we present results for W(Q m ) as a function of TO for T — 1 K, T — 15 K, and \i = —20 K for three 
values of magnetic field. Each curve shows steps as the various peaks beyond the Drude-like structure in Re tjj(£l) 
vs f2 are included in the integral defining W{VL m ) as Q TO is increased. These steps are most pronounced for the solid 
(blue) curve with B — 1.5 T. For this field value the peaks in the long dashed curve of Fig. Q] will be shifted upwards 
by a factor of ~ 1.94, so that the first drop corresponds to integrating over the first negative peak in the equivalent of 
the the long dashed (red) curve of Fig. [T] while the second drop corresponds to integrating over the second negative 
peak etc. The drops in the other curves can be similarly traced and fall at smaller values of O m because of the \f~B~ 
scaling of the peak positions in Ret#(0) vs il. To end we note that W(Cl m ) could also show upward rather than 
downward steps if we had chosen a larger value of between Mi and M 2 . In this case as an example, the first peak 
in Rei#(f2) is positive, so it adds to the low Q Drude-like contribution. 

VI. CONCLUSION 

We have studied the modification of the usual optical sum rules brought about by the particular band structure of 
graphenc namely, its compensated semimetal aspects and the Dirac nature of its quasiparticles. The existence of 
two bands within the same Brillouin zone leads to interband as well as intraband optical transitions both of which 
enter the optical sums. 

We considered both the usual /-sum rule on the real part of the longitudinal conductivity a xx (tt) and an equivalent 
sum rule recently introduced by Drew and Coleman [f| involving the Hall conductivity. From a formal point of view 
we find that care must be used in introducing the electromagnetic vector potential. It should be introduced into the 
initial Hamiltonian through a Peierls phase factor rather than after the diagonalization. The /-sum rule is shown to 
be proportional to ~ —0.39 of the kinetic energy per atom rather than the usual minus one-half value, familiar for a 
square lattice with nearest-neighbor hopping only. This difference reflects the particularity of the honeycomb lattice 
with nearest-neighbor hopping between the two distinct A and B sublattices. 

For small chemical potential \fj,\ « T. the longitudinal sum rule displays a T 3 temperature law rather than 
conventional T 2 law of free electron theory. Also for >Twe predict a \fi\ 3 dependence on chemical potential for 
the deviation from the T = /.i = reference case. The first temperature correction in this case goes as |/i|T 2 , recovering 
the T 2 law of free electron theory. The \/i\ 3 cc |V^| 3 / 2 dependence should be observable in field-effect devices. 

The sum rule for the Hall conductivity involves the real part Rei#(0) of the Hall angle (|1.3[) . Its integral (|1.4[) 
is equal to the Hall frequency ujh, and in the free electron metals ujh coincides with the cyclotron frequency. For 
graphene we find instead that ujh is given by Eq. (|4.9p . Although this equation can be formally written in the form 
(I4.10[) which involves the "relativistic" cyclotron mass m c = \n\/vp, expression (|4.10p is different from the conventional 
one. There are two modifications. The cyclotron mass m c varies as the square root of the carrier density, and uh in 
Eq. (|4.10[) is also proportional to the ratio (fi/t) 3 . 

In our previous paper [l7| we considered partial spectral weight W xx (Cl m ) given by Eq. (|5.1[) for Re <j xx (Q) and found 
stcplike structures corresponding to peaks in Rea xx (fl). Here we considered the partial spectral weight W(Q m ;B,n) 
given by Eq. (|5.2|) for Ret#(f2). Interesting stcplike structures are also found in which the steps can go up or down 
depending on the variable used to display the weight W(Q m ; B, fx). Specifically we analyzed its dependence on the 
value of the external magnetic field B and chemical potential for various fixed values of Sl m . We also considered the 
dependence of W(0, m ; B, fi) on f2 m for fixed // and B values. The rich pattern of behavior can be traced back to the 
frequency dependence of Ret# (fi) which shows a Drude-like peak at small f2 followed by a series of peaks which fall 
at some energy between the interband peaks of cr xx (fl), where both its real and imaginary parts are small. Also the 
peaks in Rei# (Q) depend on the sign and value of chemical potential. 

We hope that these specific predictions for the behavior of Re£#(f2), the corresponding partial spectral weight 
W(il m ; B, fi), and the sum rules can be verified experimentally. 
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APPENDIX A: DIAMAGNETIC TERM 

Substituting the Fourier transform (|2.1ip into Eq. (|2.7D and summing over the lattice sites, we obtain 



K 2 



op 



d 2 k 



n 2 



E 



l BZ 

d 2 k r& 



(2tt) 2 



a CT (k)V k5 6 CT (k) +h.c. 



BZ 



(2n) 2 



d 2 (f>(k) f d 2 0*(k) t v n\ 

^-Tt (k). + T CT (k) + QjTv-n (k)<x-T CT (k) 



(Al) 



Accordingly the thermal average (T a p) is easily expressed in terms of the GF (|2 . 15[) as 



2Ve 2 



n 2 



BZ 



d 2 k 

(27T) 2 



T E « 

n=— oo 
:>2, 



tr[<7 + G(iw n , k)] + tr[o-_G(iw„, k)] ^' ; ' k ' 



T -> 0. 



dk a dkp ' dk a dkp 

Calculating the trace we obtain for the diagonal component of the stress tensor (T Qj g) = (T QQ )5 Qj a that 

d 2 k 1 



(7oa)^_2e 2 T ~ 

h 2 ^ 



y V n ^J BZ {2v) 2 {uj n -in) 2 +e 2 {k)V y " ; dk 



d 2 

(t>0^)^TJ<P* : (k)+c.c 



The sum over the Matsubara frequencies converges irrespectively of infinitesimally small r, and we find 

/ " \n ,. I, (Hi - , I-, (till I , ,(Ur 

<% 2 



'" k ' [n f ( £ (k))-n F (- e (k))] ^(k)^0*(k)+c.c. 



(A2) 



V ^ y BZ (2tt) 2 2e(k) 

Writing 0(k) = e(k)e iv, ( k ' we can recast the last equation in the form ()3.5p . The derivatives in the brackets in Eq. 
are easily calculated using the explicit expression for 4>(k): 

m &f {k)+C - C - =1 2 E (<Kk)|^*(k)+c.c.) =-ye 2 (k). 



(A3) 
(A4) 



(A5) 



Substituting Eq. (| A5|) into Eq. (|A4|) we arrive at the final expression (|3.6I) . 



APPENDIX B: HALL TERM 



Here we calculate the commutator I xv — ([J x (t = 0),J y (t = 0)]) of the paramagnetic currents J x (t = 0) and 
J y (t = 0) that arc defined in Eq. (|2.18|) . Since this commutator is nonzero only in the presence of a magnetic field, 
the corresponding current density operator io(n) has to be taken for a finite vector potential A: 

ite 



Ja(n) 



(Bl) 



8,0 



where the operator 



rp _ e i«(p+e/cA)/fi_ 



(B2) 



13 



Using the anticommutation relations for a nj(T and & njCT operators we calculate the commutator 

\j x (n),j y (m)} = 6 1 * nm J2 S * S 'y ( T 6 T -6' - T &' T -s) «n, ff + bt >a (T_ S T S , - T_ S ,T S ) b n . a ) . (B3) 

" S,6',<r 

Provided that the magnetic flux per unit cell is far less than a flux quantum h/e, we may expand the operator Tg 

(B4) 



T 4 «l + i%(p a +~A a ). 



Then using the commutator 



[D a -D b 



p a + -A a ,p b + -A b 

c c 



-i — Be ab , a, b = 1,2 
c 



(B5) 



we obtain 



T&T_ S , - T S ,T_ & = ±[S a 5 b D a D b _ 5' a S b D a D b ] = 5 -A(-i)^Bz ab , 
where the sum over the dummies a, b is implied. Now the sum over 5,6' in Eq. (|B3[) can be evaluated as follows 



(B6) 



a 

E S xS' y S 7 S' K e 1K = e JK E 5 X S 7 E 5' y 8' K = e xy — 



6.6' 



where in the second equality we used the relation 



E fa)«fah 



(B7) 



(B8) 



A=l,2.3 



Substituting the commutator (|B3|) into the expression I XyV = ^ n m {[j x (t — 0, n) 7 j y (t = 0, m)]) and utilizing Eqs. (|B6 
and (|B7|) we arrive at the representation 



_ ,e 2 a 4 t 2 heB \^, ( t /.t h \\ 

11. a 



4ft 4 c 



Thus to complete the calculation of I xy we have to find the thermal average 

11. (T 

As in Appendix / is expressed in terms of the GF f|2.15p as 
„ _ „ Jbz V-^Y 



where 



Gn(iw„,k) = G 22 (iw n ,k) 



1 



1 



iuj n + [i — e(k) iuj n + fj, + e(k) 



(B9) 



(BIO) 



(Bll) 



(B12) 



are the diagonal in the sublattice space components of the GF (12.15|) . In contrast to the convergent Matsubara sums 
(|A2[) and (|A3[) in the off-diagonal components of the GF (|2.15p . regularization of the Matsubara sum (|B1 1|1 by the 
factor e~ WnT is crucial, because it would diverge otherwise. The factor e~ iuJ " T regularizes the sum over Matsubara 
frequencies differently depending on the sign of r: 



lim TV 



iuj n + fi-uj 2 L y n 2 \n F (uj)-l, r ->• +0. 



(B13) 



Normally to count the number of particles [3l| one takes the limit r — > —0. One can check, however, that such a 
prescription leads to a contradiction when calculating the integral (|B11[) , because it results in / which is not odd in 
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p, while the LHS of Eq. (|1.4p is odd in p. This contradiction is resolved if one takes into account that the first term 
of the GF (|B12[) describes electrons, so that the prescription r — * —0 has to be used, while the second term of the 
GF (|B12|) describes holes and the prescription r — > +0 is necessary. Then we obtain 

I = 2V J Bz -^[n F (e(k))+n F (-e(k)) - 1], (B14) 

which is odd in p. The quantity (|B14[) describes the carrier imbalance p — n e — rih, where n e and rih are the densities 
of electrons and holes, respectively. It was considered in the Dirac approximation in Appendix C of Ref. 20, and here 
we use two expressions (|4.7p and (|4.1ip for p. 
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FIG. 1: (Color online) Real part of the Hall angle, Retail) at B — 0.4 T, temperature T = 5K, scattering rate F = 10 K for 
four values of chemical potential. Long dashed (red) fx — — 20 K, dash-dotted (black) fx = — 300 K, solid (blue) n = — 400 K, 
short dashed (green) fi = —500 K. 
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FIG. 2: (Color online) Partial Hall spectral weight W(Q m ) as a function of field B at /i = — 20 K, temperature T = IK, 
scattering rate T — 15K for four values of Q m . Long dashed (red) f2 m = 300 cm - , dash-dotted (black) Q, m = 700 cm -1 , solid 
(blue) Q, m = 1000 cm -1 , short dashed (green) fl m = 1200 cm -1 . 
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FIG. 3: (Color online) Partial Hall spectral weight W(fl m ) as a function of /i at B = 0.4 K, temperature T = IK, scattering 
rate F = 15K for four values of Q m . Long dashed (red) Q m = 300 cm -1 , dash-dotted (black) fi m = 450 cm -1 , solid (blue) 
fl m = 700 cm -1 , short dashed (green) Q, m = 1000 cm -1 . 
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FIG. 4: (Color online) Partial Hall spectral weight W(Qm) as a function of the cutoff fi m at fi = — 20 K, T = IK, scattering 
rate F = 15K for three values of B. Long dashed (red) B = 0.5 T, dash-dotted (black) B = 0.8 T, solid (blue) B = 1.5 T. 



